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We present experimental results on eigenfunctions of a 
wave chaotic system in the continuous crossover regime be- 
tween time-reversal symmetric and time-reversal symmetry- 
broken states. The statistical properties of the eigenfunc- 
tions of a two-dimensional microwave resonator are analyzed 
as a function of an experimentally determined time-reversal 
symmetry breaking parameter. We test four theories of one- 
point eigenfunction statistics and introduce a new theory re- 
lating the one-point and two-point statistical properties in 
the crossover regime. We also find a universal correlation be- 
tween the one-point and two-point statistical parameters for 
the crossover eigenfunctions. 

PACS numbers: 05.45. Mz, 03.65. Sq, 11.30. Er, 84.40.Az 



Many complex quantum systems whose underlying 
classical behavior is chaotic can be described by treat- 
ing their Hamiltonian matrix elements as random num- 
bers which fluctuate around zero with a Gaussian distri- 
bution. There are universal statistical properties of the 
eigenvalues and eigenfunctions of these random matrices 
which depend only on the symmetries of the Hamiltonian. 
For instance, random matrix theory has been shown ti 
be consistent with the statistical properties of nuclei,! 
moleculesp and two-dimensional quantum dots.cln3 In 
the case of two-dimensional quantum mechanical systems 
with classically chaotic dynamics, there is a direct rela- 
tion between the statistics of measured conductance val- 
ues through quantum dots in the Coulomb blockade limit 
and the statistics ofjamphtudes of chaotic electron waves 
in confined systems H~l3 

When time-reversal symmetry is present, wave chaotic 
systems have statistical properties described by a Gaus=. 
sian Orthogonal Ensemble (GOE) of random matrices.Q 
As a magnetic field is applied, time-reversal symme- 
try (TRS) is lost in these systems, and the statisti- 
cal properties are described by a Gaussian Unitary En- 
semble (GUE) of random matrices. However, it is 
found in the semiclassical regime that the evolution be- 
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tween these types of symmetry is continuous and tlrrt 
a broad crossover of intermediate statistics exists.El"EJ 
It has been proposed that careful measurements of this 
crossover behavior provides-.a demanding test of the 
random matrix hypothesis J13 and we perform such a 
test in this paper. Prior experimental evidence for the 
existence of a crossover regime comes from statistical 
properties of eigenvalue spacings, which showed indica- 
tions of a progression from GOE to GUE statistics as 
a function |jjif-|tiiiie-reversal symmetry breaking (TRSB) 
parameter .Blld'O 

Here we address the evolution of eigenfunctions of 
semiclassical wave chaotic systems from the TRS to the 
TRSB limits. A considerable theoretical literature has 
developed proposing detailed descriptions of eigenvector 
statistics in the crossover regime, although little exper- 
imental data is available to test these theories. These 
theories only treat the evolution of the one-point sta- 
tistical property of eigenfunction distribution, P(|^p), 
which quantifies the degree of prohahility density, 
fluctuations in the eigenfunctions .qEjEZI No eigenfunc- 
tion imaging experiment has explicitly demonstrated the 
crossover of eigenfunction statistics from GOE to GUE 
symmetry, to our knowledge. In addition, no work has 
addressed the question of which of the theories of one- 
point eigenvector distribution function best describes the 
crossover regime and no work has addressed the crossover 
properties of the two-point correlation functions. Hence 
this work forms an important testing ground for theo- 
ries of wave chaotic systems and quantum dots based on 
random matrix theory, supersymmetry, and semiclassical 
techniques. 

The experimental arrangement used to create and 
measure the wave, chaotic eigenfunctions has been de- 
scribed previously.oE3il3 Briefly, a two-dimensional mi- 
crowave cavity with walls deflning a non-integrablc infl- 
nite square well potential is used to simulate the solu- 
tions to the two-dirpensional Schrodinger equation in the 
semi-classical limit .Ej The cavity is a symmetpjj-reduced 
bow-tie with dimensions shown in Fig. |l|(a)Jl3 A mag- 
netized microwave ferrite incorporated into the cavity 
is used to break TRS. We image the probability den- 
sity by measuring the electric fields in the 
standing wave pattern of the resonator and using the 
analogy between the Helmholtz and Schrodinger equa- 
tions in two dimensions .E3 Previous results have estab- 
lished that GOEEil and GUE statistical properties of both 
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eiff envalueJi'i and eieenfunctionslii are seen in the limit 
of zero and laege non-reciprocal phase shift in the mag- 
netized ferrite|lil respectively. 

We have found that the non-reciprocal property of the 
ferrite, hence the degree of TRSB, is a function of fre- 
quency of [the eigenmode in a relatively narrow range of 
frequency.llS This fortuitous property creates a series of 
eigenmodes of similar energy spanning the GOE to GUE 
crossover regime. Making use of this property, we iden- 
tified a crossover of A3 spectral rigidity statistics from 
the GOE to GUE limits in earlier experimental work.Eil 
In this paper we systematically examine the crossover 
eigenf unctions, and for the first time quantify the degree 
of TRSB with an experimentally determined parameter. 
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FIG. 1. Experimentally determined two-dimensional 
probability amplitude eigenfunctions plotted as logio[ 
\'^{x,y)f Acavity], derived from microwave resonator eigen- 
modes. The cavity has a symmetry-reduced bow-tie shape of 
area AcavUy and is loaded with a magnetized ferrite on the left 
side. The images are taken at a) 11.73 GHz (GOE), b) 10.79 
GHz (crossover), and c) 11.05 GHz (GUE). The antenna are 
located at (18.0, 15.5) and (45.5, 26.5) cm. 

Figure |] shows three eigenmodes |^(a;,y)p of the mi- 
crowave resonator with different degrees of TRS, but of 
similar energy. Figure |^(a) shows a GOE-limit eigen- 
mode. Fig. Mc) shows a GUE-limit eigenmode, while 



Fig. |l|(b) shows a mode of intermediate statistics. The 
probability amplitude is plotted as logiod^P^cauiii/) and 
presented in three shades to accentuate the differences 
between the GOE and GUE characteristics of the eigen- 
functions. Note that the GOE mode is distinguished by 
the tall sharp fluctuations of (dark areas), and the 
abundance of low j^fp regions (white areas) between the 
spikeso The GUE mode on the other hand has smaller 
fluctuations of l^'p and the peaks are more spread out, 
showing an aiamdance of intermediate value j^&p regions 
(gray areas) Hj These qualitative observations df^ sub- 
stantiated by more rigorous statistical analysiaiS pre- 
sented below. The crossover eigenfunction (Fig. |l|(b)) 
is a mixture of co-existing regions showing GOE- like and 
GUE-like pttpnerties. Although we see no obvious signs 
of scarringr^r^l in any of the eigenfunctions, it may make 
a contribution to the statistical measures of crossover be- 
havior presented below .CJ 




FIG. 2. Probability amplitude distribution functions plot- 



ted as P(logio(i')) vs logio(t;), where v = 1*1 Acavity Solid 
lines represent theoretical distribution functions in the GOE, 
GUE, and crossover regimes described by the Zyczkowski and 
Lenz theory Pb{v). Also shown are averages of distribution 
functions for data in the crossover regime. 

As mentioned above, the distribution of j^'p values is 
a simple means of identifying the symmetry of the cor- 
responding Hamiltonian. In the GOE limit one finds 
the Porter-Ihomas distribution function of eigenfunc- 
tion values,E3 Pgoe{v) — (l/\/27rt;)e^"/^, where v ~ 
Acavity , and Acavity IS the area of the cavity. In the 
GUE li«i±-.the distribution function is simply Pgue{v) 
— e^".EZlc3 To enhance the subtle differences between 
these two distribution fupctions, we plot P(logio(w)) 
versus logio(w) in Fig. |^.EJ Note that P gue {^ogio{v)) 
peaks at a higher value and falls off more quickly, while 
PGO-E(logio(i^)) has a shorter and broader distribution. 
From this it is clear that GOE eigenfunctions have more 
large (black in Fig. |^) and small (white) j^'p fluctua- 
tions, consistent with our qualitative observations of Fig. 
|l|. Also shown in Fig. ^ are distribution functions of 
eigenmodes in the crossover regime. Each data set is an 
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average of 5 eigenfunctions of similar degree of TRSB 
(as discussed below). There is a smooth variation of 
the distdbution of eigenfunction fluctuations as TRS is 
brokenE3 The exact form of this variation is not predicted 
by theory. 

We have found four theories for the GOE 
GUE crossover probability density distribution function 
Pcrossoveriy^ that duc to Brickmann, et al. (Pp^iA, 
P e (1,2)),B Zyczkowski and Lenz (P^), b e (1,2)),EI 
Sommers and lida (Pe(w), £■£ (0, cx))),llZl and Fal'ko and 
Efetov (Px(w), X G (0, oo)).a The subscript denotes the 
symmetry-breaking parameter used by the authors. In 
each case, the lower limit of the TRSB parameter (/3, 6, 
£, X) corresponds to GOE statistics, and the upper limit 
to GUE statistics. All of the theories agree in these two 
limiting cases, but they disagree foiiintermediate param- 
eter values in the crossover regime.EZI For instance Ppiv) 
= (w), a generalized chi-square distribution suggested 
by Brickmann, et shows a maximum of the distribu- 
tion function Pp{v) at u = 1 for all values of /3, whereas 
the other theories show a peak in Pb,e,x{v) which occurs 
for u < 1 in the crossover regime (see e.g. the solid lines 
in Fig. I for Pbiv)). 

It should be noted that none of the theories for the 
crossover behavior is precisely aplicable to this exper- 
iment. The theory of Brickmann, et al. proposed an 
interpolation of the chi-square distribution with no mi- 
croscopic justification. Zyczkowski and Lenz performed a 
more detailed calculation in which two chi-square distri- 
butions were convolved together because it was assumed 
that the real and imaginary parts of the random matrix 
elements were independently fluctuating. Sommers and 
lida addressed the issue of properly treating the fluctua- 
tions of the matrix elements, but had to perform an en- 
ergy averaging to arrive at the final distribution function. 
The work of Fal'ko and Efetov uses a model appropriate 
for disordered systems and was not explicitly derived for 
ballistic billiard systems such as ours. (However the re- 
sults here suggest that the theory is also applicable to 
such systems.) Despite the fact that none of these the- 
ories exactly corresponds to the experimental situation, 
they are remarkably accurate in their description of the 
data. 

The probability distribution functions P{v) for 64 ex- 
perimental eigenmodes have been fit to the four theories 
of crossover statistics. It is found that the Zyczkowski 
and Lenz theory produces the best fits for these eigen- 
functions, although the difference in fit quality for the 
other theories is not statistically significant. Figure ^ 
shows averaged distribution functions for groups of 5 
eigenfunctions with similar values of the Zyczkowski and 
Lenz crossover parameter 6, as well as best fits to the 
averaged Pb{v). Notice that significant changes in P{v) 
already occur for a small deviation of b from 1 , consistent 
with Fal'ko and Efetov's prediction that a small amount 
of magnetic flux quickly moves the system away from 
GOE statistics^ Also note that the peak of the experi- 



mental distribution function occurs at w < 1 for crossover 
eigenfunctions, demonstrating that the generalized x| (v) 
distribution is not a correct description of the crossover 
data. |— I 

These theories, as well as random matrix theory,ll3 do 
not predict how the crossover parameters evolves with 
non-reciprocal phase shift in the ferrite, and do not even 
relate their crossover parameters to those of other theo- 
ries. Therefore we proceed empirically and define a sim- 
ple experimental measure of the degree of TRSB. We 
have noticed a strong asymmetry of the forward (Si2(/)) 
and reverse (S2i(/)) complex transmission coefficients of 
the microwave cavity for TRSB eigenmodesjij and cor- 
related this asymmetry with the degree of TRSB derived 
from statistical analysis of the associated eigenfunctions. 
We define the experimental time-reversal asymmetry 
parameter A = J i\\Si2\ - \ S 2i\\) df / J (|5i2| + |^2i|)rf/, 
where the integrals are carried out over one resonant 
peak (between neighboring minima in |S'i2(/)|) in the 
frequency domain. This parameter is easily jeyaluated, 
does not require an image of the eigenmodej^ and (as 
shown below) can be considered a measure of the "mag- 
netic fiux" through the cavity causing TRSB. 



2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.0 



0.8 



□ P 

• b 

A (2E+1)/(e+1) 

X (2X+1)/(X+1 



0.0 0.1 0.2 0.3 

Asymmetry Parameter, A 



0.4 



FIG. 3. Plot of crossover parameters /3, b, (2e + 1) /(e + 1), 
and {2X + 1)/{X + 1) versus the experimentally determined 
TRSB parameter A. The GOE limit is in the lower left, while 
the GUE limit is in the upper right of the figure. The error 
bars show the typical standard deviation of the data points. 

We find that the one-point correlation function pa- 
rameters (/3, b, e, X) are strongly correlated with our 
experimental measure of TRSB, A. Figure ^ shows 
the four distribution function crossover parameters /3, 6, 
{2e + l)/(e + 1), {2X + 1)/(X + 1) plotted vs. A, for 
groups of similar eigenfunctions. Note that the param- 
eters £ and X have been mapped to the interval (1,2) 
in an ad hoc way, simply for comparison with the other 
statistical parameters. We see that all four statistical pa- 
rameters describe a smooth and universal transition from 
GOE to GUE statistics as the asymmetry of the trans- 
mission characteristic. A, increases.Eil Note that there is 
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initially a linear increase of the crossover parameters with 
A, suggesting that they all give a measure of the "mag- 
netic flux" producing the GOE GUE crossover. The 
reduction of the parameter values beyond A — 0.25 is 
likely due to difficulties in properly calculating A for the 
highly distorted microwave resonance curves Si2{f) en- 
countered in the GUE limit. 

One can also employ two-point statistical correlation 
functions to quantify the crossover from TRS to TRSB 
behavior of the eigenfunctions. It has bepj shown that 
C(fcr) = (|^'(0)«'(r)|2) - 1 + cJ^{kr)jB where k is 
the wavenumber of the eigenmode, Jo is the zero-order 
Bessel functioru-jawd c = 2 for GOE and c = 1 for GUE 
eigenfunctions. ta-Ea Eigenfunctions with a smaller value 
of c show a more "smeared out" appearance (Fig. |l|(c)). 
We find that the value of c changes smoothly between 2 
and 1 for the experimental crossover eigenfunctions. 

To investigate the self-consistency of these results, we 
have examined the relationship between the one-point 
correlation function parameters (/3, 6, (2e -t- l)/(e + 1), 
{2X -I- 1)/(A^ -I- 1)) and the two-point parameter c (de- 
termined by fitting C(fcr)), as shown in Fig. ^. Again 
we notice that there is a universal behavior of the o: 
point crossover parameters in their dependence on c. 
This demonstrates that the crossover behavior is shared 
by both single point and two-point statistical properties 
of the eigenfunctions, in agreement with prediction 
of the supersymmctric nonlinear a model.c 



2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.0 
0.8 



ak|e GOE corresponds to 
k + 6-k zero (to make ^ real). 



1 



0.6 




A (2e+1)/(e+1 
X (2X+1)/(X+1 



0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 

Two-Point Correlation Parameter, C 



2.4 



FIG. 4. Plot of crossover parameters /3, b, (2e -I- l)/(e + 1), 
and {2X + 1)/{X + 1) versus the two-point correlation func- 
tion parameter c. The GUE limit is in the upper left, while 
the GOE limit is in the lower right of the figure. The error 
bars show the typical standard deviation of the data points. 
The inset shows a correlation between the one-point statisti- 
cal parameter b, and the two-point parameter c, in the GOE 
to GUE crossover regime. 

To develop a relationship between one-point and two- 
point eigenfunotipn properties, consider a superposition 
of plane waves,E3 ^' ~ o-ke^^^ , with fixed wavenum- 
ber, k — |k|, but with random directions and ampli- 



tudes. Writing Ok = 
|ak| = |a-k| and Ak = 
while GUE corresponds to Ak random in [0,27r]. Al 
lowing general Ak in the crossover regime leads to the 
two-point result for C{kr) given above with c = 1 -I- 
I ^ |ak||a_k| exp(iAfe)/^ lokPP, and to the one point 
result of Zyczkowski and Lenz, P6(v), with c = 1-1- 

(t^)^- '^'^ ^^^^ ^^^^ prediction, we plot c — 1 versus 
(t^)^ as an inset to Fig. ||. The data points fit to a 
straight line of slope close to 1, demonstrating remark- 
able agreement between the one-point and two-point sta- 
tistical parameters in the crossover regime. 

To summarize, we have established both GOE (TRS) 
and GUE (TRSB) properties of wave chaotic eigenfunc- 
tions and shown that there is a continuum of eigenfunc- 
tions with intermediate statistics between the two time- 
reversal symmetry states. We have found that three 
of the four theories of one-point eigenfunction statis- 
tics P(f ) in the crossover regime adequately describe our 
data. We have introduced a simple experimental quan- 
tity which measures the degree of TRSB, and shown 
that the statistical properties of the eigenfunctions evolve 
smoothly with an increase of this parameter. Finally we 
demonstrate that the one-point and two-point correlation 
functions describe the smooth crossover in a consistent 
and universal manner, and that a single TRSB parameter 
must describe the crossover behavior. 
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